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ABSTRACT 

During three balloon flights of a 1 sr ionization- 
chamber/Cerenkov counter detector system, we have measured 
the atmospheric attenuation, flux, and charge composition 
of cosmic-ray nuclei with 16 < Z < 30 and rigidity greater 
than 4.5 GV. 

The attenuation mean-f ree-path in air of VH (20 < Z < 30) 

2 

nuclei is found to be 19.7 ± 1.6 g/cm , a value somewhat 

greater than the best previous measurement. The attenuation 

2 

mean-f ree-path of iron is found to be 15.6 ± 2.2 g/cm , 
consistent with predictions of geometric cross-section 
formulae. 

We measure an absolute flux of VH nuclei 10 to 20% higher 
than earlier experiments at similar geomagnetic cutoff and 
level of solar activity. The relative abundances of even- 
charged nuclei are found to be in good agreement with results 
of other recent high-resolution counter experiments. 

We calculate that our observed cosmic ray chemical 
composition implies relative abundances at the cosmic-ray 
source of Ca/Fe = 0.12 + 0.04 and S/Fe = 0.14 ± 0.05. The 
results are consistent with all other elements of charge 
between 16 and 26 being absent at the source and being 
produced by cosmic-ray fragmentation in interstellar hydrogen. 
The results show a cosmic-ray source deficient relative to 
the solar system in both sulfur and argon. 
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1. INTRODUCTION 

We have measured the flux and composition of cosmic-ray 
nuclei with atomic number (Z) greater than 15 and rigidity 
greater than 4.5 GV. These measurements were made in a series 
of three balloon flights during 1970 with a large-area detector 
using ionization chambers and a Cerenkov counter. In this 
paper we present results on nuclei with 16 < Z < 30. Our 
measurements of VVH (Z > 30) nuclei will be presented in a 
later publication. 

o 

With our very large geometrical factor (0.92 m ster) , 
we have observed more than 25,000 VH (20 < Z < 30) nuclei, 
several times the sum of all previous investigations. Al- 
though our charge resolution is not as good as that of some 
other investigators (Webber e_t cH_. , 1971; Casse ejb aK , 1971) 
we have used our excellent statistics to measure properties 
not easily accessible to these other experiments, including 
the atmospheric attenuation of these nuclei, and their 
absolute flux. In addition we have measured relative 
abundances of even-Z nuclei and extrapolated these abundances 
back to the cosmic ray source. Our calculations demonstrate 
that either a simple slab model or a steady-state exponential- 
path-length model of cosmic-ray propagation can fit the 
observed abundances of nuclei with 16 < Z < 30 with nearly 
the same source abundances. 
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Preliminary results of these three flights and of a later 
flight at lower geomagnetic cutoff have been presented at the 
Twelfth International Conference on Cosmic Rays (Binns et al . , 
1971). Further detail on this experiment and the data analysis 
may be found in Mewaldt, (1971), and a more detailed detector 
system description appears in Epstein et al . , (1971). 
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2. THE DETECTOR SYSTEM AND BALLOON FLIGHTS 

A schematic diagram of the detector system is shown in 
Figure 1. The principal components are the two pulse-type 
ionization counters and the Lucite Cerenkov counter. The two 
paral 1 el -pi ate ionization counters are 1 m in area, 5.1 cm 
thick, and filled with a mixture of 90% argon and 10% methane 
at a pressure of one atmosphere. The Cerenkov radiator is 
a 1 m sheet of 0.6-cm thick ultraviolet-transmitting (UVT), 
Lucite, with light collection by total internal reflection. 
Each of the four edges of the radiator is viewed by a 12.7-cm 
photomultiplier (RCA 4525) through an adiabatic light pipe 
constructed of UVT Lucite. The outputs of the four phototubes 
are summed. Approximately 2.5 photoelectrons result from the 
perpendicular traversal of a relativistic single-charged 
particle. 

The particle trajectories are given by a hodoscope con- 
sisting of two crossed layers of scintillator strips above 
the other counters, and two crossed layers below. Each layer 
consists of ten strips, giving 9.5 cm resolution. A six-fold 
coincidence among the four hodoscope layers and the two ion 
chambers triggers the recording of an event. The ion chamber 
discriminators determine the experiment threshold and are 
typically set at a signal level corresponding to a vertically 
incident relativistic nucleus of charge 16. For each event. 
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the signals from each of the ion chambers and from the Cerenkov 
counter are pulse height analysed, and the trajectory of the 
particle is indicated by recording which hodoscope-stri p dis- 
criminators were triggered. 

There are four guard counters, one on each side of the 
detector square. The guard counters are not in active anti- 
coincidence, but the outputs of the four guard counter dis- 
criminators are recorded for each particle triggering the 
system . 

This paper treats data gathered in three high altitude 
balloon flights of this detector system (Table I). The residual 
atmosphere, geomagnetic cutoff, and neutron monitor values 
averaged over the entire float are presented here for comparison 
however, in applying corrections to the data, hourly averages 


were used. 



BALLOON FLIGHTS 
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3. CORRECTIONS TO RAW DATA 


The raw data were modified by correcting recorded pulse 
heights for temporal and spatial variations in detector 
response and by eliminating the small fraction of ambiguous 
events . 

Three corrections were applied to the ionization chamber 
pulse heights. Corrections for variations of the temperature 
and pressure of the gas were made using pulse height spectra 
from the internal Americium - 241 al pha -par ti cl e calibration 
source. Hodoscope data indicating position in the chambers 
of each particle's trajectory were used to correct for areal 
non-uniformities in the chamber response (typically 1 to 2 %). 
Hodoscope data indicating zenith angle (e) of each particle's 
trajectory permitted correction of each pulse height for the 
sec 6 dependence of the path length in the gas. 

The Cerenkov counter pulse heights were corrected for 
non -uni f ormi ti es in light collection which were determined 
empirically using cosmic-ray iron nuclei. The zenith angle 
dependence of Cerenkov pulse height was much stronger than 
sec 0, in general agreement with the performance of a similar 
counter described by Ormes and Bal asubrahmanyan (1970). More 
significantly, the variation of light collection with velocity 
is considerably greater than the variation of light production. 
In addition, there are systematic azimuthal variations. We 
can explain the measured response of the counter in terms of 
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the amount of light produced, the fraction that is internally 
reflected, and a semiempirically -determi ned light collection 
efficiency (Mewaldt, 1971). 

In selecting events for further analysis, we have elimi- 
nated only those events for which charge identification was 
highly ambiguous. Our criteria rejected 4% of all recorded 
events. Appropriate small corrections were made in determining 
absolute fluxes. In eliminating the ambiguous events we 
avoided criteria which might preferenti al ly eliminate some 
charges. The acceptance criteria were: 

1. Ion Chamber Agreement --Even ts in which the two ion 
chambers disagreed by more than three charge units were re- 
jected. The number of such events (~1%) is consistent with 
the expected number of nuclear interactions in the material 
between the chambers . 

2. Ion Chamber-Cerenko v Cons i s tenc.y --Events were re- 
jected if the ratio of Cerenkov to ionization signals was 
more than three standard deviations greater than the maximum 
theoretical value. These were less than 0.2% of all events. 

3. Multiple Hodoscope Events --In about 7.5% of all 
events, two adjacent hodoscope strips were triggered in at 
least one hodoscope layer, in many cases the result of knock- 
on electrons. These events were normally accepted. However, 
when two non-adjacent strips were triggered, (-2.5% of all 
events) the trajectory was ambiguous, and the events were 

re jec ted . 
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4. Guard Events --0ne or more guard counters was trig- 
gered in less than 2% of the events. Since these events 
appear to be normal in other respects, they were rejected 
only if two or more guards fired, or a simultaneous multiple 
hodoscope occurred, amounting to less than 0.5% of all events. 

We note that the 4% event rejection is significantly 
less than the fraction of events that undergo nuclear inter- 
actions within the detector system. Indeed, many of these 
interactions are not distinguishable from normal events. 

At the energies of the particles detected in these flights, 
interaction secondaries closely preserve the direction of 
the primary nucleus, and our relatively coarse hodoscope 
resolution does not resolve the individual fragments. To 
correct for interactions in the material above the ion 
chambers (1.2 g/cm ), we extrapolated the observed spectrum 
through an equivalent thickness of atmosphere. Interactions 
below the ion chambers have almost no effect on the charge 
assignment, as it is determined almost solely by the ion 
chamber measurements. 
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4. PARTICLE IDENTIFICATION 

Both the rate of energy-loss of a particle (measured 

by the ionization chamber) and its Cerenkov light product ion 

2 

are proportional to Z , but the two quantities have different 
functional dependence on the particle's velocity, so it is 
possible to use these measurements to determine both the 
particle's charge and velocity. Figure 2 shows a cross-plot 
of the square root of the ionization (/T) vs the square root 
of the Cerenkov light intensity (/C) for even-Z nuclei. We 
have included the measured response of our ionization counters 
in the relativistic rise region (Fernandez, 1972) and the 
response of our Cerenkov counter at 0 = 30°. For the balloon 
flights reported here, the events are expected to fall to the 
right of the 4.5 GV geomagnetic cutoff indicated on the cross- 
plot. Note that charge and energy identification are ambiguous 
in the relativistic rise region where, for example, an iron 
nucleus of 20 GeV/nucleon produces the same I and C as a 
2.0 GeV/nucleon nickel nucleus. Because of the steepness of 
the cosmic ray energy spectrum, this ambiguity is a problem, 
only when the lower charge is much more abundant than the 
higher charge. 

While it is apparent that for high cutoff flights, the 
Cerenkov counter provides only limited improvement over charge 
identification based on ionization alone, the Cerenkov counter 
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is very useful in insuring that rare particles of very high 
ionization are actually relativistic VVH nuclei, and not slow 
VH nuclei . 

Figure 3 shows a /T vs /C cross-plot for about 2000 

events with 20 < 6 < 50°. The predicted response of the 

counters agrees well with the data. Notice that there are 

virtually no events lying significantly to the right (high 

C/I ratio) of the predicted response. There are, however , 

a substantial number of events (~6%) that appear to have 

energies less than that allowed by the geomagnetic cutoff 

(low C/I ratio). These events originate from particles that 

2 

undergo nuclear interactions in the 0.7 g/cm of the Lucite 
Cerenkov counter after passing through the ion chambers. 

Because ionization is nearly independent of energy for 
rigidity greater than 4.5 GV, a good first approximation to 
charge identification is given by Z a /T. An improvement 
in charge resolution of about 20% was achieved for events 
due to particles with energy between 1 and 2 GeV/nucleon 
where the ratio C/I gives an unambigous, charge-independent 
measure of energy. However, for events with low Cerenkov 
signals, and those in the relativistic rise region (T > 2 GeV/ 
nucleon), a velocity determination becomes impossible, and 
charge identification was based on ionization alone, using 
the mean of the two ion chambers. 
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Figure 4 shows a charge spectrum obtained from -10,000 
events with 0 < 30°, collected at an atmospheric depth of 
4.5 g/cm and a vertical geomagnetic cutoff near 4.5 GV . 

The iron peak resolution (~1.5 charge units fwhm) is typical 
of the charge resolution obtained in this experiment. 

The charge scale in Fig. 4, derived by assuming that 
the prominent peak is centered at iron, agrees to within 
0.1 charge units with a calculated scale based on the ioniza- 
tion chamber response to the in-flight alpha-particle calibra- 
tion source (Epstein ejt aK , 1971). This agreement confirms 
the generally accepted assumption that the cosmic ray flux 
has a significant peak at Z = 26. 

In order to identify the relative contributions of the 
individual elements to this charge spectrum, a least squares 
fit was performed. The fit used a calculated pulse-height 
distribution for any individual element which took into 
account the measured response of our ion chambers in the 
relativistic-rise region and left the charge resolution as 
a free parameter. The resulting fit gave unambiguous results 
for the relative abundances of elements of even charge. The 
elements of odd charge are less abundant, and our resolution 
did not permit determination of their abundances. Nickel 
and Zinc (Z = 28 and 30), though of even charge, were obscured 
by the relativistic rise of the much more abundant iron, so 
their abundances were not well determined. 
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5. ATTENUATION AND FRAGMENTATION IN THE ATMOSPHERE 

Taking advantage' of our large geometrical factor, we 
have measured the atmospheric attenuation of VH nuclei using 
data from the ascent portions of the balloon flights. 

Charge spectra from each flight were obtained in each 

2 

of seven depth intervals over the range 4.0 < x < 22.5 g/cm . 
The thickness of atmosphere penetrated by each incident 
particle was determined from its trajectory in the detector 
and the altitude of the balloon at the time of arrival; thus 
particles penetrating any given thickness of atmosphere were 
observed from various zenith angles at various balloon alti- 
tudes. The effective area x solid-angle x time factor for 
each depth interval was calculated by numerical integration 
of the detector's differential geometry factor over the 
balloon ascent trajectory. The data from the three flights 

were normal i zed to a 4.5 GV vertical geomagnetic cutoff, as- 
- 1 4 

suming a P * integral rigidity spectrum, and to a Mt. 
Washington neutron monitor rate of 2128 using the helium 
regression curve of Webber (1967). 

Figure 5 shows the depth dependences of the Fe and of 
the VH fluxes along with the results of linear least squares 
fits. The data points have been plotted at the mean atmos- 
pheric depth traversed within each depth interval, as computed 
from the experimental data. Effects due to variation of 
geomagnetic cutoff with incident direction are negligible 
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because each data point represents an average over all azimuths 

and a wide range of zenith angles. 

The fit to the data yields an attenuation mean free path 

2 

for iron of Ap e = 15.6 ± 2.2 g/cm , in good agreement with 

the interaction mean free paths computed from the semi -empi ri cal 

o 

formulae of Daniel and Durgaprasad (1962), 14.2 g/cm ; Bradt 

and Peters (1950), 14.0 g/cm^; and Cleghorn et al . (1968), 

2 2 
15.5 g/cm . Our value for (19.7 ± 1.8 g/cm ) is notice- 

2 

ably greater than the value of 16.5 ± 2 g/cm obtained by 
Webber and Ormes (1967), although the uncertainties overlap. 

The value of Ay H is greater than Ap e primarily because 
of the non-zero probability of VH nucleus remaining in the 
VH group after an interaction (p VH VH > 0). Neglecting small 
differences in the geometrical cross sections of iron nuclei 
and VH nuclei, we calculate 

PVH.VH 1 - °- 21 * °' 13 (,) 

This is consistent with the value 0.26 ± 0.05 obtained by 
Waddi ngton ( 1 960) . 

Extrapolation to the top of the atmosphere of fluxes 
of nuclei with Z < 26 cannot assume a simple exponential 
attenuation, because these lighter nuclei are produced 
copiously in the atmosphere from iron fragmentation. It 
is convenient, instead, to extrapolate the ratio R(x) = 
N 2 (x)/N-|(x) where is the flux of nuclei in group i and 
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x is the atmospheric depth. If group 1 (e.g. iron) has 
negligible contribution from atmospheric fragmentation of 
heavier nuclei and group 2 consists of lighter nuclei whose 
only significant fragmentation source is group 1, it can be 
shown (Mewaldt, 1971) that R(x) varies linearly with x for 
x << a i ( A 2 — A i ) as 

P ijAo Ao-Ai 

R(x) . R(o> ♦ tR(o) ♦ (2) 

where a ^ is the attenuation mean free path and p^ is the 
probability that an interacting particle of group 1 which 
is thrown out of group 1 by the interaction falls into 
group 2. 

The observed depth dependence of each appropriate flux 
ratio (Fig. 6) was fit according to the above equation. The 
resulting slopes were used to extrapolate flux ratios observed 
at float altitude to the top of the atmosphere. (Results are 
given in the next section.) 

An estimate of the atmospheric fragmentation spectrum of 
the iron group (25 < Z < 30) was obtained by investigating the 
depth dependence of the ratio N ( 2 5 < Z < 30)/N( Z Q < Z < 24) 
for 18 < Z Q < 24, using Eq. (2). Although the accuracy of 
this determination was limited by statistical uncertainties, 
the data was consistent with a flat fragmentation spectrum, 
with p 25 3 q 20 24 - 0.25 ± 0.12 for the fragmentation of 
25 < Z < 30 nuclei into 20 < Z < 24 nuclei. 
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6. CHARGE COMPOSITION 

For each flight, a least squares fit to the data was 
used to determine the charge spectrum at the detector.. 

This spectrum was then extrapolated to the top of the 
atmosphere. In Table II we summarize our results on the 
charge composition of nuclei with 16 < Z < 30, where the 
results at the top of the atmosphere represent an average 
over the two flights at P c = 4.5 GV. For comparison, we 
show also the charge composition at detector level (4.1 g/cm ) 
from flight 1970 - 2 . The quoted uncertainties include sta- 
tistical and charge identification uncertainties and un- 
certainties in the extrapolation to the top of the atmosphere. 

Except for the ratios involving sulfur (Z = 16), these 
results were derived using only events with e < 30° and the 
data were extrapolated to the top of the atmosphere as 
described in the previous section. Because of the level 
of the iron chamber discriminators, we observed nuclei with 
Z < 17 only at wide angles. The MH/VH ratio [ ( 1 6- 1 9 ) / ( 20-30 ) ] 
was derived using data from flight 1970-2 in the angular range 
40 < 9 < 60°, and extrapolated to the top of the atmosphere 
(~6 g/cm ) using = 22.5 g/cm and p^ H ^ = 0.15 (Freier 
and Waddington, 1968). The sulfur abundance was extrapolated 
to the top of the atmosphere using the fragmentation parameters 
of Long ( 1 968) . 
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TABLE II 

RELATIVE ABUNDANCES 


p 

Charge Ratio Top of Atmosphere 4.1 g/cm 

Average of two Flights Flight 1970-2 


24/26 

0.15 

+ 

0.03 

0.16 

± 

0.02 

22/26 

0.14 

+ 

0.03 

0.16 

± 

0.02 

20/26 

0.21 

+ 

0.04 

0.22 

± 

0.03 

18/26 

0.13 

± 

0.03 

0.14 

± 

0.02 

16/26 

0.27 

± 

0.06 + 

0.33 

± 

0.06* 

( 21 -24) /( 25-30) 

0.32 

± 

0.04 

0.35 

± 

0.02 

( 20-24 )/( 25-30) 

0.45 

± 

0.05 

0.50 

± 

0.03 

( 1 8-20) / (21 -30) 

0.24 

± 

0.03 

0.26 

± 

0.02 

(1 6-19)/( 20-30) 

0.31 

± 

0.05 f 

0.36 

± 

0.05* 


^'Flight 1 970-2 only 
*At ~6 g/cm^ 
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Figure 7 compares our results for even-Z nuclei with 
those of other recent experiments of good charge resolution. 
Our results agree well with those of Webber e_t aj_. (1971) 
and Casse et al . (1971). 
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7. ABSOLUTE FLUX 

The absolute flux of relativistic VH nuclei has been 
determined by a number of workers at various geomagnetic 
cutoffs and levels of solar activity. The results with the 
greatest statistical weight until this experiment are those 
of Webber and Ormes (1967), Freier and Waddington (1968), 
and von Rosenvinge ejt aK ( 1 968, see von Rosenvinge, 1 970 , 
for a more complete summary of this work). In general, our 
fluxes are consistently higher than those obtained by these 
other workers. 

Figure 8 and Table III summarize our results. In Table 
III we have selected data from our flights obtained at the 
same vertical geomagnetic cutoff as the flights of these 
other workers. All measurements have been extrapolated to 

2 

the top of the atmosphere using the value a vh = 19.7+1.8 g/cm 
determined in this experiment. Normalization to a Mt. Washing- 
ton neutron monitor rate of 2350 was accomplished using the 
helium regression curve of Webber (1967). The uncertainties 
are those given by the authors. In the case of von Rosenvinge 
(1970), they include only statistical uncertainties. Errors 
quoted on our results include effects of uncertainties in 
charge identification, in corrections for the few rejected 
events, and in extrapolation to the top of the atmosphere, 
as well as the relatively small statistical errors. We note 
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the general sel f-consi stency of our measurements over three 
flights. We also obtain consistent VH flux measurements 
over the area of our counter, and at various zenith angles, 
taking into account the appropriate atmospheric attenuation. 
This degree of consistency allows upper limits of only a 
few percent to be placed on the number of out of geometry 
events that may have been accidentally included in the 
analysis. We have applied a number of small corrections 
for the effects of interactions and losses within our 
counter, but the product of these amounts to less than 4%, 
within our quoted uncertainty. 

We note from Figure 8 that our VH fluxes are system- 
atically about 20% higher than those of von Rosenvinge. 

Our disagreement with Freier and Waddington is smaller. 
Although the fluxes of these other experiments are derived 
from only a few hundred particles, the disagreement is 
appparently not of statistical origin. 

As the VH charge composition derived from the Freier 
and Waddington 4.9 GV flight (see Waddington ejt al_. , 1970) 
is significantly different from that found by other workers, 
it is conceivable that charge identification problems have 
caused them to find a low VH flux. 

Two of the flights of von Rosenvinge involved systematic 
corrections of the order of 100% for VH nuclei, the result 
of losses of events within the detector system (von Rosen- 
vinge, 1970). As no estimate of the uncertainty in this 
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correction is included, it is difficult to assess the serious- 
ness of the disparity. A large error in the determination of 
the residual atmosphere is apparently ruled out by the fact 
that the disagreement between von Rosenvinge and ourselves 
spans three flights with each detector. As the flight of 
von Rosenvinge at 4.1 GV occurred at a neutron monitor rate 
(2125) equal to that of our flight 1970-4 (2128), we can 
apparently eliminate the possibility of different levels 
of solar modulation. 

It is of interest to look at these discrepancies in 
terms of the charge composition of the VH group. Von 
Rosenvinge obtains a ( 20- 24 ) / ( 25-30 ) ratio of 0.61 at the 
top of the atmosphere from two flights at 4.35 and 4.9 GV. 

We obtain at 4.5 GV a ( 20-24 )/( 25-30) ratio of 0.45 ± 0.05, 
in agreement with the recent measurements of Webber et al . 
(1971) and Casse et_ a_l_. (1971). We find that our measurements 
of the (20-24) group flux are in general agreement with von 
Rosenvinge, thus the discrepancy between our VH flux measure- 
ment and that of von Rosenvinge appears to be concentrated 
in the iron group. Our value for the MH(16 < Z < 19) flux 
(0.14 ± 0.03 m~ sr sec" ). is also in apparent agreement 
with von Rosenvinge's value. 
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8. COMPOSITION AT SOURCE 


In this section we present cosmic ray source composi- 
tions derived from our measurements. We have calculated 
the composition which would be observed at earth assuming 
various source compositions and propagation parameters. 

The experimentally observed spectrum was then used to place 
restrictions on the source composition and mode of propagation. 

In the "leaky box" model (Goeckler and Jokipii, 1969), 
it is assumed that the cosmic ray sources are uniformly 
distributed in space and time, such that an equilibrium 
is established between creation and loss. The particles 
are assumed to be confined to the galaxy, although there 
is a finite probability for escape which can be characterized 
by a mean life against escape, T, commonly referred to as 
the "leakage lifetime". This steady state model is mathe- 
matically equivalent to an exponential path length distribution 
( Dav i s , 1 960 ) . 

In the steady state the number density, N., of particles 

J 

of type j with a uniform continuous source Q., is then described 

J 

by 


w .i (x) 
A .• 


"i ( *> 


l 

i >J 


N i (x) 


i J 


+ Q 


0 


where the leakage length, X, is pcT, with p the interstellar 
gas d e n s i ty . 
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The fragmentation parameters (p. .) were calculated from 

' J 

the semiempi ri cal cross-section formula of Shapiro and Silber- 
berg (1970) for a representative energy of 2.5 GeV/n, and the 
A for a nucleus of atomic number A from 




A = Mp/a 

where 

% 

= proton mass 

and 

a 

■ .(-.A 1 ' 3 )* 

wi th 

r o 

= 1.2 x 10' 1 3 cm 


We have assumed a pure hydrogen medium. 

Our neglect of energy loss processes is justified pro- 
vided ( dE/dx ) /E << j- + ^ (Gloeckler and Jokipii, 1969). 

This inequality is satisfied for the relativistic particles 
we observe. (For Fe 3 ^ nuclei at P c = 4.5 GV, (dE/dx)/E= 0.04 
and + j = 0.6). The error due to neglect of energy-loss 
is further minimized since we treat only charge ratios of 
i ntegral fluxes, over a limited charge region. 

In Figure 9 we compare the detailed predictions of the 
steady state calculation with our observations of even-Z 
nuclei. We have adjusted the leakage length so that the 
observed events with 21 < Z < 24 are accounted for by inter- 
stellar fragmentation ,of iron. (It is reasonable to assume 
a very low source abundance for these elements in view of 
their very low abundance in the solar system, 0.02 of iron. 
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(Cameron, 1968)) Assuming a pure iron source (solid line) 
large discrepancies appear at Z = 16 and Z = 20, implying 
a substantial source abundance of these nuclei. This con- 
clusion is not surprising in view of the pronounced peaks 
in the experimental charge spectrum (Fig. 4) at calcium 
and sulfur. 

To explain the observed abundances, we require a calcium 
to iron ratio at the source of 0.12 ± 0.04. For sulfur we 
find ( s / Fe ) S0urce = 0.14 ± 0.05. The dotted line in Figure 
10 shows the effect of including this amount of calcium and 
sulfur at the source. Notice that fragmentation of the 
additional calcium also improves the agreement at argon 
(Z = 18) . 

2 

The leakage length X = 4 g/cm implied by this compari- 
son is consistent with values that are found to fit the L/M 
ratio and VVH nuclei. (Gloeckler and Jokipii, 1969; Cowsik 
ejtaJ_.,1970; Shapi ro e_t aJL , 1970; Fowler ejt a_K , 1970). 

We have also performed propagation calculations for the 

simple "slab model". The observed (21 -24)/(25-30) ratio sets 

2 

an upper limit of X Q = 1.85 ± 0.2 g/cm to the amount of 
interstellar hydrogen traversed. Although this small value 
of X Q rules out the slab as a generally applicable model 
(because the observed ratio of light to medium nuclei requires 
X Q of 3 to 5 g/cm ) it is still instructive to derive source 
abundances over a 1 imi ted i nterval of charges using this model. 
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Table IV summarizes the inferred source composition 
and compares the result with the Cameron (1968) solar system 
abundances. These cosmic ray results depend only weakly on 
the (21 -24)/(25-30) ratio at the source. For example, if 
(21 -24) / ( 25-30) source = 0.1 we find Ca/Fe = 0.15 ± 0.04 and 
S/Fe * 0.18 ± 0.05 on the steady state model. We conclude 
that for either the slab (X Q = 2 g/cm ) or the steady state 
model (X = 4 g/cm ) the observed high energy charge spectrum 
for 16 < Z < 26 is consistent with a source composed primarily 
of iron, with smaller amounts of calcium and sulfur. In fact 
the relative abundances of sulfur, calcium, and iron at the 
source are essentially the same in either model. The remaining 
abundances are consistent with a fragmentation origin and 
these elements are not required in any large concentration 
at the source, although abundances of a few percent that of 
iron are not excluded. 

Several other authors have recently considered the 
source abundances of cosmic ray nuclei in this charge region, 
including Shapiro and Silberberg (1970), von Rosenvinge (1970), 
Webber et_ £l_. (1971), Cartwright e_t aj_. (1971), and Shapiro 
et al . (1971). There now appears to be general agreement 
among all workers on the abundance of sulfur at the source. 
Shapiro e^t aj_. (1971), using a weighted average of data 
from many experiments published prior to 1971, have recently 
revised their calculations (Shapiro and Silberberg, (1970), 
and now include small and roughly equal amounts of argon and 
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TABLE IV 

ABUNDANCES RELATIVE TO IRON 


Cosmic Ray Source 
(slab model , 

X = 1.85 g/cm 2 ) 

Cosmic Ray Source 
( 1 eaky box model , 
X - 4 g/cm 2 ) 

Solar System 
(Cameron, 1968) 


Ca 1 c i urn 


0.13 ± 0.04 


0.12 ± 0.04 


0 .08 


Sulfur 
0.16 ± 0.05 

0.14 ± 0.05 
0.57 
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calcium at the source, (A/Fe = Ca/Fe = 0.09 ± 0.09). The 
results of our investigation indicate a slightly higher 
source abundance of calcium, and we definitely find A < Ca , 
in general agreement with von Rosenvinge (1970), and Webber 
e_t aK (1971). Earlier experimental measurements showed a 
sizeable Cr/Fe ratio at both high and low energies (Mathiesen 
et al . , 1968, Haddington e_t aj_. , 1 970 ; Lezniak ejt al . , 1970; 
Dayton ejt aj_. , 1 970; Price ejt aj_. , 1 970 ; Garcia-Munoz and 
Simpson, 1970), forcing the conclusion that substantial 
chromium existed at the source (Shapiro and Silberberg, 

1970, von Rosenvinge, 1970, Longerfelter and Ramaty, 1971). 

More recent measurements with better resolution and statistics 
show a much lower flux of chromium (Webber ejt a_K , 1971; 

Casse e_t aj_. , 1971). The chromium abundance observed in 
this experiment is low enough to conclude that it is princi- 
pally of secondarry origin, in agreement with the conclusions 
of Shapiro et al . (1971), and Webber et al . (1971). At low 
energies, the recent satellite measurements of Cartwright 
et al . (1971) at 20-300 MeV/n imply a source composition 
with Ca/Fe = 0.1 and S/Fe = 0.16 in agreement with that 
found here. 

We conclude that the trend of the most recent experimental 
and theoretical work favors finite source abundances of calcium 
and sulfur, with all other elements (with Z > 16) less than 
10% of Fe at the source. Future work will concentrate on 
some of these less abundant elements. 
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Comparison between our cosmic-ray source results and the 

Cameron solar system abundances indicates a slight cosmic-ray 

excess in calcium and a marked cosmic-ray deficiency in argon 

and sulfur. We note, however, that the solar system argon 

abundance is very uncertain (as with other volatile elements) 

and Cameron's abundance for A is simply an interpolation 

32 40 

between observed abundances of S and Ca . The low cosmic- 
ray sulfur abundance appears to represent a real disparity 
with the solar system. 

We consider it unlikely that the differences we find 
between the cosmic-ray source and solar-system composition 
are a result of propagation effects. More likely, the ex- 
planation is to be found in the environmental conditions 
under which these nuclei were synthesized. It now appears 
that solar-system elements in this charge region were syn- 
thesized under explosive conditions in massive stars, as 
evidenced by the dramatic success of explosive nucleosynthesis 
theory in matching the solar-system abundances (Arnett and 
Clayton, 1970). With'the cosmic ray elemental (although not 
yet isotopic) abundances now known to an accuracy comparable 
to the solar-system abundances, it should soon be possible 
to identify the environment where these nuclei originated. 
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FIGURE CAPTIONS 

1. Schematic cross-section of the detector system. 

2. Square root of ionization vs square root of Cerenkov 
pulse-heights for various nuclei. 

3. Cross-plot of square root of ionization vs square-root 
of Cerenkov pulse-height for 2000 events. 

4. Charge histogram. At charge 16 and below the number 
of events is depressed due to electronic threshold. 

Arrows above charge 30 indicate single events. 

5. Atmospheric attenuation of VH nuclei and of iron nuclei. 
Symbols indicate flights. a-1 970-2 , X-1970-3, 0-1970-4. 

6. Atmospheric depth dependence of various charge ratios. 

The data points are averaged over the three flights. 

7. Comparison of abundance measurements for even-Z nuclei 
at the top of the atmosphere. 

8. Comparison of VH flux measurements. All measurements 

have been extrapolated to the top of the atmosphere 

2 

using Ay^ = 19.7 g/cm and normalized to a Mt. Washington 
neutron monitor rate of 2350. 

9. Comparison between observed abundances and abundances 
calculated in steady state model for a pure iron source 
and for a source with calcium and sulfur abundances 12% 
and 14% of iron respectively. 
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